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Abstract 

Buckminsterfullerene (Ceo) was recently confirmed to be the largest 
molecule identified in space. However, it remains unclear how, and where 
this molecule is formed. It is generally believed that Ceo is formed from 
the build up of small carbonaceous compounds, in the hot and dense en- 
velopes of evolved stars. Analyzing infrared observations, obtained by 
Spitzer and Herschel, we found that Ceo is efficiently formed in the ten- 
uous and cold environment of an interstellar cloud illuminated by strong 
ultraviolet (UV) radiation fields. This implies that another formation 
pathway, efficient at low densities, must exist. Based on recent laboratory 
and theoretical studies, we argue that Polycyclic Aromatic Elydrocarbons 
are converted into graphene, and subsequently Ceo, under UV irradiation 
from massive stars. This shows that alternative - top-down - routes are 
key to understanding the organic inventory in space. 



1 Introduction 

The mid-infrared spectra of a variety of astrophysical objects are dominated by 
band emission (strongest at 3.3, 6.2, 7.7, 8.6 and 11.2 fim) attributed to car- 
bonaceous macromolecules, i.e. Polycyclic Aromatic Hydrocarbons [T] (PAHs). 
These molecules are large (30-100 C atoms), abundant (~ 5% of the elemental 
carbon), and their ionization plays a key role in the energy balance of gas in 
the interstellar medium and in protoplanetary disk. In addition to PAH bands, 
infrared signatures observed at 7.0, 8.5, 17.4 and 19.0 /im have been reported 



recently [U |3] , and found to coincide precisely with the emission of buckmin- 
sterfullerene (Ceo) [1], a cage-like carbon molecule. This detection heralds the 
presence of a rich organic inventory and chemistry in space. However, observed 
abundances of Cgo challenge the standard ion-molecule or grain-surface chem- 
istry formation routes, which build up molecules from small to large in the 
interstellar medium. For that reason, it has been suggested that Cgo is formed 
in the hot and dense envelopes of evolved stars O EJ [7j in processes alike to 
those found in sooty environments O El [TOl [H] , and eventually, is ejected in 
space. Yet, this scenario faces the problem that it has a limited efficiency [6]. 
PAHs and Cgo are known to co-exist in the interstellar medium [3] , however, so 
far, the connection between PAHs and Cgo -and in particular the possibility to 
go from one compound to the other in space- has not been investigated. In this 
paper, we present a study of PAH and Cgo chemical evolution in the NGC 7023 
nebula, using Spitzer |12j and Herschel 13 infrared observations. 

1.1 Observational results 

1.2 Infrared observations of the NGC 7023 nebula 

Earlier Spitzer observations of the NGC 7023 reflection nebula have revealed 
a chemical evolution of PAHs: deep in the cloud, emission is dominated by 
PAH clusters, which evaporate into free-flying PAHs when exposed to the UV 
radiation from the star [Ml [151 [E] ■ There, gaseous PAHs are, in turn, ionized. 
While the neutral PAHs are dominated by zig-zag edges - as demonstrated by 
the strong C-H solo out-of-plane modes -, the ions have an arm-chair molecular 
structure, characterized by strong duo out-of-plane modes [T7]. In regions closest 
to the star, the presence of Ceo in the neutral state is evidenced by Spitzer 
observations (Fig. 1). New Herschel observations provide a measurement of 
dust emission in the same region, at high angular resolution (Fig. 1). This 
measurement can be used to derive the integrated intensity radiated by the 
nebula which can be used as a calibrator, to convert the Spitzer observations 
of the PAHs and Ceo bands into absolute chemical abundances of these species, 
allowing a quantitive study of PAH and Ceo chemical evolution. 

1.3 Measurement of the far infrared integrated intensity 
of dust emission in the nebula 

The far infrared integrated intensity Ifir was extracted by fitting the spectral 
energy distribution (SED) at each position in the cross cut shown in Fig. 1. For 
these positions we have used the brightnesses as measured by Herschel PACS 
Poglitsch et al. (20 10)1 and SPIRE [22] photometers. We have used the 70 and 
160 jxm channels of PACS and the 250 /im channel of SPIRE. This data is 
presented in details in [2S|. The modified blackbody function fit to these SEDs 
containing 3 spectral points is defined by: 

I{X,T) = K/\P xB{\T), (1) 
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where K is a. scaling parameter, A is the spectral index, and B{X, T) is the 
planck function with A the wavelength and T the temperature. We have used a 
constant value of 1.8 for /3 so that only T and K are free parameters. The results 
of the fits to the observations are shown in Figure SI. The temperatures we 
have derived from the fit of the data range between 25 and 30 K. These values 
are in agreement with those of [2 3) for the same region. The peaking position 
of the modified blackbody function moves to shorter wavelengths (i.e. the grain 
temperature increases) when getting closer to the star, implying that we are 
indeed tracing matter inside the cavity and not material behind on the line of 
sight. The far infrared integrated intensity, IfiRi is then derived by integrating 
I{X,T) over frequencies. 

1.4 Measurement of the integrated intensity of PAH and 
Cgo emission in the nebula 

The integrated intensity of PAH emission, Ipah, is measured by fitting a PAH 
emission model to the observed Spitzer mid infrared low resolution spectrum. 
This data covers the 5 to 14 fj,m. range where most of the emission occurs. 
Because this spectral range contains the emission due to the vibration of both 
C-C and C-H bonds, it is insensitive to ionization of PAHs. An example of this 
fitting model is presented in [21]. This tool provides the integrated intensity 
in the PAH bands as an output and takes care of continuum subtraction and 
extinction correction. The Cgo integrated intensity (Iceo) extraction needs to 
be done separately for each band. We first measure the integrated intensity 
in the 19 and 17.4 band. The first step consists in extracting the intensity in 
the 19.0 ^m band (/ig.o) by fitting a gaussian and subtracting a local linear 
continuum. The 17.4 /xm band is contaminated by PAH emission, but this 
can be removed effectively. The contribution of PAHs to the 17.4 band can be 
estimated in this way: in the outer regions of the nebula, the 17.4 /xm band 
is 100% due to PAHs (since no Cgo is detected there) and so is the 16.4 ^m 
band. The 16.4 and 17.4 /im band of PAHs are known to correlate. Indeed, in 
the outer regions of the nebula we find ly/^^^ ~ /i6.4 x 0.35. Therefore, over 
the whole nebula, we can estimate the intensity of the 17.4 /im band due to 
Ceo by Ii'^°4 ~ Ii7A — Ii6a x 0.35. Figure S2 shows the result of this process 
on the map. The 7.0 fim band is faint and usually hard to detect. In the 
regions closest to the star where both the 19 and 7.0 features are observed, 
we can calibrate the ratio of /ig.o/^T.o- It is found to be relatively stable (at 
least in this small zone close to the star) and of the order of 0.4. Therefore 
we use /7.0 — 0.4 x /19.0. The 8.5 fim band is undetectable because of the 
strong PAH band present at 8.6 /zm, so we use the ratio provided in |3] for 5 eV 
photons that is /g 5 = 0.4 x /ig.o- The total Cgo IR. emission is hence given by 
Ic^o = ^19-0 + hjA - h&A X 0.35 + /19.0 X 0.4 + /19.0 x 0.4. 
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1.5 Derivation of PAH and Ceo abundances 



The method to derive the abundance of carbon locked in PAHs and Cgo from 
I PAH, Iceo ^^'^ I FIR is presented in details in |Tielens(2005)| . The fraction 
of carbon locked in PAHs and Ceo (respectively Jq^" and fQ^°) per atom of 
interstellar hydrogen are given by: 



and 



fS- = <,.n.C-^)^J^. (3) 



a- 



uv 



1 - (Rc'eo)' 



where 



RpAH^ ^ , (4) 

^PAH + ^Ceo ' ^FIR 



and 



IpAH + ICeo + IPIR 

'^uv'^ and cr^^" are the UV absorption cross sections of PAHs and Cgo- Fol- 
lowing [Tielens(2005")] , we adopt tif^^^ = 7 x 10"^^ cm^ per C atom. There are 



no detailed measurements of the UV absorption cross section of Cgo, following 
[Sj we adopt the same value as for PAHs. 



1.6 Evidence of Cgo formation in NGC 7023 

The results of the abundance variations within the nebula as derived using the 
above method are shown in Fig. 2. The Cgo abundance in the nebula is seen to 
increase from 1.4 x 10"''% to 1.5 x 10~^% of the elemental carbon abundance 
when approaching the star (Fig. 2). On the other hand, the abundance of 
Polycyclic Aromatic Hydrocarbons (PAHs) is seen to decrease, from 7.0% to 
1.8% of the carbon (Fig. 2). This shows that Cgo is being formed in the nebula 
while PAHs are being destroyed or processed. The correlation of these variations 
with the increasing UV field strongly suggests that UV photons control Cgo 
formation and PAH processing and destruction. 



2 The top-down model 

2.1 Proposed scenario for the formation of Ceo in the ISM 

The formation of Cgo in the interstellar medium is unexpected and has several 
implications in our understanding of the formation process of this molecule. In 
the laboratory, the main formation route invoked for Ceo is the build-up from 
atomic carbon, small carbon clusters or rings [HI [HI 113 IH] ■ In space, such pro- 
cesses are efficient in the hot (1500 K) and dense {nn > 10^^ cm~^) envelopes 
of evolved stars 1^. In NGC 7023, the gas is ~ 7 orders of magnitude less dense 
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(see [23] and SI text), making this aggregation processes inefficient. Instead, we 
invoke photochemical processing of PAHs as an important route to form Ceo- 
Upon UV irradiation, several channels for fragmentation can be open depend- 
ing on excitation energy; e.g., H-loss and C2H2 loss ^7\. However, experiments 
on small PAHs have shown that H-loss is by far the dominant channel [3S] 
and that complete dehydrogenation i.e. graphene [22] formation, is then the 
outcome of the UV photolysis process [IHl Tielens(2005)[ 1501 IM] . In a second 



step, carbon loss followed by pentagon formation initiates the curling of the 
molecule [32]. We envision that this is followed by migration of the pentagons 
within the molecule, leading to the zipping-up of the open edges forming the 
closed fullerene [T^ (Fig- 3 and SI Movie). Graphene formation through PAH 
photolysis can give rise to a rich chemistry. Besides the route towards fullerenes 
outlined above, fragmentation towards small cages, rings, and chains may also 
result. The relative importance of isomerization and fragmentation will deter- 
mine the carbon inventory delivered by the photo-chemical evolution of PAHs 
in space (Fig. 3). 

2.2 PAH and Graphene stability in space 

Guided by experimental [TS] and theoretical [50] studies we have evaluated the 
stability of PAHs against dehydrogenation, and Cgo formation, in conditions 
appropriate for NGC 7023. 

2.2.1 Graphene formation 

Schematically, the fragmentation process can be written as, 

PAH + hv ^ PAH* PAH_H + H , (6) 

where PAH* is the excited species which can stabilize through emission of IR 
photons or through fragmentation and PAH_^f is a dehydrogenated PAH rad- 
ical. There are various ways to evaluate the unimolecular dissociation rate 
constant for this process (cf., |Tielens(2005)] sect 6.4). We wih follow [531 and 
write the rate constant in Arrhenius form, 

k (E) = K {Te) exp [-EjkT,] , (7) 

where Tg is an effective excitation temperature, Eq-, the Arrhenius energy de- 
scribing the process, and the preexponential factor, kg, depends on the interac- 
tion potential (in the reverse reaction). For PAHs, the internal excitation of the 
vibrational modes after the absorption of a UV photon can be well described 

by, 



fE{eY) 



OA 



T™ = 2000 l^-^ j , (8) 

where is the number of C-atoms and E is internal energy in eV ( [Tielens(2005)| 
p. 184). Because typically the energy involved in these reactions is a fair fraction 
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of the total energy in the system, a correction has to be made to this excitation 
temperature. The finite heat bath correction results in Tielens(2005)[ 155] . 



E 

The preexponential factor can be set equal to, 



= T„, ( 1 - 0.2 § ) . (9) 



(10) 



ko = —r- exp 
n 



' ■ R 



with AS the entropy change for which we will adopt 5 cal/K [5D]. ko is then 
typically, ~ 3 x 10^^ s~^. The Arrhenius energy parameter, Eo, cannot be 
easily evaluated from theoretical calculations ([34])- Here, we use a fit to the 
experimental fragmentation studies on small PAHs (< 24 C-atoms; [28|), which 
resuhs in Eo is 3.3 eV ( [Ticlcns(20'05)] p. 204). The probability for dissoci- 
ation depends then on the competition between fragmentation and IR photon 
emission, 

where kir{E) is the IR emission rate at an internal energy E. For a highly 
excited PAH, kir{E) is ~ 1 s~^. The total fragmentation rate is then, 

Pd {E) fc™ [E) , (12) 

where kuv (E) is the absorption rate of UV photons with energy, E. 

The photochemically driven H-loss is balanced by reactions of atomic hydro- 
gen with dehydrogenated PAHs; viz., 

PAH_H + H ^ PAH + hiy . (13) 

The rate of this reaction has been measured to be ka = 1-5 x 10~^° cm^ 
for a number of small PAHs [301 ISS] ■ We can define the dissociation parameter, 
ip — kaUH /kfrag- With kuv = 7 X 10~^° A^c Go (Go is the flux of UV photons 
in units of the interstellar Habing field [35], we have ij} ~ 0.2nH /NcGoPd{E). 
The hydrogen coverage of interstellar PAHs is a very sensitive function of tp and 
a small increase of ip can change PAHs from fully hydrogenated to graphene 
[Tielens(2005)[ [55] . This is illustrated for the circumovalene, C66H20, in Figure 
S3. Thus, PAHs are fully hydrogenated if tp is much less than 1 and fully 
dehydrogenated if ip is much larger than 1. Hence, 

Go 0.2 ^ kUE)\ ^^^^ 



UH No V k{E) 

provide a critical relation for the transformation of PAHs into graphene. We 
have evaluated this relation assuming an absorbed UV photon energy of 10 eV 
(Fig. 4). 
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2.2.2 Fullerene formation 



After formation of graphene, UV photoabsorption can lead to loss of carbon from 
the skeleton. This fragmentation process competes with stabilization through 
IR photon emission. The reaction rate is again given by eqn. [7j We adopt AS 
is 5 cal/K, resulting in fco ~ 3 x 10^^^ s^^, but the exact value is not critical. 
For Eo, we have adopted the Arrhenius energy (3.65 eV [37]) derived from 
experiments on the C-loss from small catacondensed PAHs [28| • The calculated 
cohesive energy of carbon in graphene is much larger, 7.4 eV/C-atom [38] but 
that refers to typically carbon inside the skeleton and carbon at the edge will 
be less strongly bound. [T^ find a theoretical value of 5.4 eV for the C atoms at 
a zig-zag edge. Moreover, theoretical cohesive energies are not a good measure 
for the Arrhenius energy in unimolecular dissociation experiments [34j . Given 
the low abundance of gas phase carbon, the reverse reaction is unimportant and 
the chemical lifetime (t^^^ = kfrag{E j) has to be compared to the dynamical 
lifetime of the region. We have evaluated this chemical lifetime for a typical 
internal energy of lOeV as a function of the size of the graphene sheet, Nc, and 
the results are shown in Fig. 4. 

2.3 Application to NGC 7023 

To evaluate graphene formation in NGC 7023 we look at PAH stability. The 
physical conditions in the nebula, characterized by parameters Gq and njj, can 
be obtained (SI Text). As shown in Fig. 4, H loss is very efficient in NGC 
7023 for PAHs up to 70 C atoms. We then evaluate the timescale for the loss of 
C atoms by the graphene flakes and compare this to the dynamical age of the 
nebula for a ~ 70 C atom PAH (Fig. 4). It appears that if we adopt Eg = 3.65 
eV, C-loss is rapid compared to the age of the nebula. 

3 Discussion 

Our models reveal that 70 C atoms PAH become unstable to graphene forma- 
tion in the studied region of NGC 7023. Much larger PAHs will survive closer to 
the star, while smaller PAHs will rapidly be completely dehydrogenated further 
away but - because of size - this will not lead to the formation of fullerenes. 
As shown in laboratory experiments [THj, graphene sheets larger than about 
70 C atoms can be transformed into fullerene, but in space this is driven by 
UV photons rather than energetic electrons. We surmise that fullerene forma- 
tion is initiated by single C atom loss [ID] at the zig-zag edge of the graphene 
flake. Quantum chemical calculations indeed show that the armchair structure 
of graphene clusters stabilizes through in plane 7r-bond formation, while the 
open-shell orbitals associated with the dangling bonds are located at the zig- 
zag edge [35 , making these atoms more labile [19] . A detailed study of graphene 
edge structure ^39], based on experimental results [ID], has indeed shown that 
single C atom loss leads to the formation of reconstructed edges bearing pen- 
tagons. These defects dramatically stress graphene flakes, inducing significant 
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positive curvature in its topology [35] . Eventually, the barrier- less migration of 
the pentagons within the hexagonal network will lead to the zipping-up of the 
flake into Ceo ([IS]) Fig- 3)- The above mechanism is severely limited by the 
physical conditions, controlling dehydrogenation, and the size of the precursors, 
which need to be above 70 C atoms. As shown by our observations, only 
about 1% of the available PAHs are converted into Ceo by this process (Fig. 2). 
The activation of this "top-down" chemistry requires the high UV fields avail- 
able near massive stars. As a corollary, high abundance of Ceo in the diffuse 
interstellar medium |41j reflects processing by massive stars, which are capable 
of dehydrogenating PAHs in the relevant size range (60-70 C atoms). PAHs 
much larger than than 70 C atoms will survive even in close encounters with 
massive stars. Smaller size PAHs ( < 50 C atoms) cannot reach the fuUerene is- 
land of stability, hence their photo-processing (Fig. 2) must be a source various 
carbon nano-compounds e.g. rings, chains, cages, bowls, tubes etc ( Fig. 3). 
Laboratory experiments have indeed shown that processing of small graphene 
constrictions can lead to carbon rings and chains [31]. Finally, we note that 
recent studies suggest that PAH molecules do not provide a satisfactory expla- 
nation to the diffuse interstellar bands [331 13S] ■ The compounds described in 
this top-down chemistry may be relevant to understand these features of the 
UV extinction curve, in the Milky- Way and in other galaxies [?5] . 

4 Conclusion 

Analyzing infrared observations of the NGC 7023 nebula we have found evidence 
that Ceo is being formed in the ISM. Classical bottom-up formation routes fail 
to explain these observations so we propose a new chemical route in which Ceo is 
formed directly by the photochemical processing of large PAH molecules. Other 
carbonaceous compounds (cages, tubes, bowls) are also the product of this pho- 
toprocessing. This route must be relevant in the interstellar medium, but may 
also be important in the inner regions of protoplanetary disks around solar-type 
stars where accretion on the young star generates intense UV fields and PAHs 
are known to be present. There, they can be converted, on timescale of '-^ 1 
Myr, into fullerenes and hydrocarbon fragments. This source of organic com- 
pounds remains to be considered in models studying the organic photochemistry 
of the interstellar medium and regions of terrestrial planet formation. It is clear 
that such studies may benefit greatly from the progress achieved in the field of 
graphene stability. 
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Figure 1: Overview of the NGC 7023 nebula. (A) Multi wavelength color coded 
view of the nebula in the infrared. Red is the emission at 70 /im observed with 
the Photodetector Array Camera end Spectrometer (PACS) onboard Herschel. 
This emission is dominated by dust. Green is the Spitzer-lnira,Ted Array Camera 
(IRAC) 8 ^m emission tracing the PAH C-C mode and blue is the IRAC 3.6 /Ltm 
emission, tracing the PAH C-H mode and stellar emission. The position of the 
intermediate mass young star HD 200775 illuminating the nebula is indicated. 
(B) Color coded image of the spatial distribution of different compounds in 
NGC 7023 : Red is the emission of dust observed with Herschel-P ACS, green 
shows the emission integrated in the 6.2 /im C-C band of PAHs observed with 
/Spzteer-Infrared Spectrograph (IRS), and blue is the Cgo emission observed with 
Spitzer-IKS integrated in the 19.0 /j,m band. The white rectangle shows the 
region on which the extraction of the PAH and Ceo abundance (Figure 2) was 
performed. (C) Spitzer IRS mid-infrared spectra taken at positions (1) (upper 
panel) and (2) (lower panel) in (B). The distance from the star at positions (1) 
and (2) are respectively ~ 35" and ^ 15". The bands of PAHs and Ceo are 
labeled in the spectra. 
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Figure 2: Abundances of Ceo (A) and PAHs (B) in NGC 7023 as a function of 
distance from the star in the cut shown in Figure 1. The red curves give the 1 
sigma uncertainty, obtained from the propagation of instrumental uncertainty 
on the extraction of integrated intensities of PAHs and Cgo bands. 
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Figure 3: Schematic representation of top-down interstellar carbon chemistry. 
(A) The chemical evolution of PAHs in the interstellar medium under the influ- 
ence of UV photons combines the effects of dehydrogenation and fragmentation 
with those of isomerization. Fully hydrogenated PAHs - injected by stars into 
the ISM - are at the top right side. Near bright stars, UV photolysis will pref- 
erentially lead to complete H-loss (e.g., the "weakest link") and the formation 
of graphene. Further fragmentation may lead to the formation of flats, rings, 
and chains. However, this process competes with isomerization to various types 
of stable intermediaries such as cages and fuUerenes. (B) Schematic illustra- 
tion of conversion of graphene into Cgo in 13 steps. Dehydrogenated PAHs, 
i.e. graphene sheets, loose carbon atoms under UV irradiation, giving rise to 
pentagonal defects (represented with orange C atoms) at the edges of the sheet. 
These defects in the hexagonal network induce curvature of the sheet, the mi- 
gration of the pentagons allows the molecule to close. Image courtesy of L. 
Cadars (www.laurecadars.com). 
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Figure 4: Photochemistry of PAHs in NGC 7023 (A) The red curve shows the 
evolution of the critical value when PAHs have lost half of their H atoms, as 
a function of the physical conditions (Go/nn where Gq is the radiation field 
in Habing units and uh is the density of H atoms in the gas in cm~'^, see 
Online Supporting Text for details) and the number of carbon atoms in the 
PAH molecule Nc- Above this line PAHs are stable against dehydrogenation. 
Below this curve, PAHs rapidly loose all their H atoms. The value of Gq/uh 
for NGC 7023 is shown, dark blue represents the "narrow" range and clear blue 
the "broad" range of values for this parameter (Online Supporting Text). (B) 
Time scale for the loss of a carbon atom by UV photon absorption as a function 
of graphene sheet size. 
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Supporting Text 



Physical conditions in NGC 7023 
General properties 

NGC 7023 is a widely studied photodisociation region. We briefly review the 

information on physical conditions available in the literature. The nebiila has a 
hourglass-shaped, low-density cavity, which was opened in the dense molecular 
cloud by the winds of the young Be star (Puente et al. 1998). The formation of 
Ceo is seen to occur inside this cavity (See Fig. 1 in the main text) in regions 
particularly close to the star (10" to 40"). In the next sections we discuss the 
values we adopt for radiation field and gas density in the cavity. 

Radiation field 

The intensity of the radiation field in the cavity. Go, can be derived based on 
the spectral type and flux of the star diluted by the square of the distance from 
the star. Doing this exercise, Joblin et al. (2010) find that the radiation field 
is Go = 2600 at a distance of 40" from the star (similar to the value found by 
Rogers et al. 1995 and Chokshi et al. 1988). Using this value and spherical 
dilution we derive the range of radiation fields covered in the cut where we study 
Ceo and PAH evolution (Figs. 1 and 2). We find Go = 3.0 x 10* at the peak 
of Geo abundance (12" from the star) and Go = 3.4 x 10^ at the minimum of 
Ceo abundance (35" from the star). Fuente et al. (1999) have reported the 
highest value of radiation field in the PDR with Go = 1 x lO"*^ at the position 
of the H2 filaments at 46" from the star. This would correspond to values of 
Go = 1.7 X 10" and Go = 1.5 x 10^ at the edges of the cut (resp. 35" and 12" 
from the star). We therefore adopt average values Go = 1.0 ± 0.7 x 10* and 
Go = 1.0 ± 0.7 X 10^ at 35 and 12" from the star from the star. 

Density in the cavity 

The density of the atomic gas in this region is difficult to derive directly, but 
can be constrained from observations of the surrounding molecular cloud. More 
specifically, the density derived from CO lines by Gerin et al. (1998) in the 
back wall of the cavity, point to a value of nf^°^ ^ 3000 cm~^. The atomic gas 
well within the cavity, is expected to be at least order of magnitude hotter, and 
hence an order of magnitude less dense if we consider pressure equilibrium, i.e. 
ng"" - 300. Rogers et al. (1995) find that "'/ng'^" = 10 - 35, which then 
imphes ng''^ = 85 - 300. Joblin et al. (2010) quote a value of 100 cm'^ for 
ng"". A more direct estimation of the column density of warm atomic gas can 
be derived from the dust emission. Since we have seen that dust temperature 
increases when getting closer to the star (see first section of this document) , this 
implies that this emission indeed comes form the cavity and not from the wall 
behind. Using the DUSTEM (Compiegne et al. 2011) model, we can reproduce 
the observed emission in the PACS 70 /xm filter at a distance of 35" from the 
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star 1 Jy/pixel or 2300 MJy.sr"-'^), for a radiation field of Go = 10^, and 
leaving the column density as a free parameter. This yields Nh = 2 x 10^^ 
cm~^ for the cavity, which with a physical size of 5.7 x lO^^cm (2x35" at 400 
pc) corresponds to n^""" = 45cm~^, somewhat lower than other estimates based 
on PDR or molecular tracers. Note that this determination does not depend 
significantly on the accuracy of determination of radiation field in the considered 
range of Go (see Compicgnc ct al 2011 Fig. 7). We keep n^°^ = 50 as lower 
limit and ng'*" = 250 as upper limit and adopt ng"" = 150 ± 100 cm"^. 

Gq/uh 

Using the adopted average density and radiation field values we can obtain 
Go/riH needed to estimate the dehydrogenation efficiency (see section "Graphene 
formation" in the manuscript). We find Gq/uh = 65 ± 45 at 35" from the star 
and Gq/uh = 650 ±450 at 12" from the star. From these numbers we define the 
"broad" and "narrow" domains of values of Gq/hh expected in the cut in NGC 
7023, respectively [20-1100] and [105-200]. These are represented graphically in 
Fig. 4 over the dehydrogenation stability curve. 



18 



Supplementeiry videos 



https : // sites.google.com/site / olivierberne / c60- formation- from-graphene 

S- Video 1: Conversion of graphene into Ceo- This video shows schematically how 
graphene can be converted into Ceo. We start with a sheet of graphene where 
carbon atoms are arranged in a hexagonal network. Under UV irradiation, 
C atoms are lost in the hexagons situated at the edge of the graphene flake 
and converted into pentagons (shown in orange). The formation of pentagons 
stresses the molecule forcing it to cin-vc. When several pentagons have been 
formed they can migrate inside the molecule. Finally, when 12 pentagons have 
been formed and when each carbon atom belongs to both a pentagon and a 
hexagon, the molecule closes into Ceo- 
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Supplementary figures 
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Figure S-1: Fitted modified black bodies to the SEDs observed with Herschel 
in the cross cut from position (1) to (2) show in Fig. 1. 
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Figure S-2: Subtraction of the contamination of PAHs in the 17.4 /xm band. 
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Figure S-3 The relative abundance of circumovalene and its partial dehydro- 
genated derivatives, CeeHn, for two values of the dissociation parameter, ip. 
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